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A  PROCEDURE  FOR  CALCULATING  THE  DEVELOPMENT  OF  TURBULENT 
BOUNDARY  LAYERS  UNDER  THE  INFLUENCE 
OF  ADVERSE  PRESSURE  GRADIENTS 


By  Kennedy  F.  Rubert  and  Jerome  Persh 


SUMMARY 


In  an  attempt  to  better  predict  turbulent-boundary-layer  develop¬ 
ment  in  aircraft  diffusers,  a  procedure  based  on  the  kinetic-energy 
equation  and  an  extended  form  of  the  momentum  equation  has  been  devised 
for  calculating  the  development  of  turbulent  boundary  layers,  in  adverse 
pressure  gradients,  for  that  class  of  flows  for  which  the  fluid  density 
at  all  points  and  total  pressure  outside  the  boundary  layer  are 
invariant.  In  the  development  of  this  procedure  an  effort  was  made  not 
only  to  arrive  at  an  analytical  form  that  would  allow  examination  of 
the  significance  of  the  physical  quantities  involved  but  also  to  achieve 
a  high  degree  of  consistency  with  the  more  recent  results  of  turbulent- 
boundary-layer  research. 

Correlations,  which  are  essential  to  the  execution  of  the  method, 
of  several  elements  of  the  basic  equations  in  terms  of  quantities 
regarded  as  controlling  parameters  are  presented.  These  correlations  - 
are  strictly  interim  empirical  relationships  both  as  to  numerical  values 
and  as  to  the  variables  involved.  It  is  to  be  expected  that,  with  the 
acquisition  of  more  data  covering  a  broader  range  of  conditions,  improve¬ 
ments  will  be  made  in  both  the  functional  nature  and  the  accuracy  of  the 
correlations.  A  comparable  and  associated  refinement  of  the  basic 
equations  is  likewise  anticipated  as  more  data  are  obtained. 

Predictions,  by  this  method,  of  turbulent-boundary-layer  develop¬ 
ment  are  compared  with  experimental  results  from  several  sources  for  a 
number  of  cases  of  flow  over  flat  plates  and  airfoils  and  in  conical 
diffusers.  In  the  range  of  boundary-layer  flow  short  of  separation,  the 
agreement  of  calculated  values  with  experimental  values  was,  in  most 
cases,  quite  satisfactory;  in  some  instances,  however,  definite  disagree¬ 
ment  between  the  calculated  and  experimental  results  was  noted.  It  is 
believed,  however,  that  good  agreement  has  been  obtained  in  enough 
instances  to  justify  continuation  of  effort  along  the  present  lines, 
particularly  with  respect  to  improvement  of  the  correlations  and  refine¬ 
ment  of  the  equations. 


2 


NACA  TN  21*78 


INTRODUCTION 


Because  of  the  dominant  part  played  by  the  boundary-layer  growth 
in  determining  diffuser  performance,,  theoretical  diffuser  analysis  is 
effectively  a  problem  in  the  calculation  of  the  boundary-layer  develop¬ 
ment.  For  aircraft  application  the  boundary  layer  is  usually  turbulent 
and  of  a  thickness  less  than  the  half-span  of  the  channel.  A  semi- 
empirical  procedure  for  the  calculation  of  the  development  of  such 
boundary  layers  in  two-dimensional  flows  has  been  given  by  Von  Doenhoff 
and  Tetervin  (reference  1).  Attempts  at  direct  application  of  this 
and  other  existing  methods  to  conical-diffuser  flows,  however,  have  led 
to  difficulties  both  in  the  accuracy  obtained  under  some  conditions  and 
in  the  physical  interpretation  of  the  basic  equations.  The  present 
research  was  undertaken  with  the  object  of  evolving  a  procedure  which 
would  facilitate  interpretation  of  the  physical  quantities  and  would  be 
applicable  to  conical-diffuser  flows.  This  present  paper  is  primarily 
concerned  with  the  progress  to  date  in  achieving  this  objective.  It 
deals  with  a  procedure  developed  for  that  class  of  flows  for  which  the 
fluid  density  at  all  points  and  total  pressure  outside  the  boundary  layer 
are  invariant,  presents  interim  empirical  correlations  of  certain 
physical  quantities  essential  to  the  execution  of  the  procedure,  and 
includes  an  examination  of  the  results  obtained  to  date. 


SYMBOLS 


D 

H 

q 

R0 


IX 


u» 


U 


X 


duct  diameter 

boundary-layer  shape  parameter  (6*/e) 
dynamic  pressure  outside  boundary  layer 


boundary-layer  Reynolds  number  based  on  momentum 
thickness  (0Up/p) 

velocity  parallel  to  surface  at  a  perpendicular  distance 
y  from  wall 

root-mean-square  fluctuating  velocity  in  axial  direction 

velocity  parallel  to  surface  at  a  distance  greater  than 
or  equal  to  6  from  wall 

longitudinal  distance  parallel  to  wall  measured  from 
initial  station 
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% 

Y 


* 


mi 


% 


P 

'To/2'! 

(V2q)s.Y. 

(to/2ci)l.t. 

^r/2^ 

ts/2(1 

Subscripts : 

exp 

corr 


duct-station  locations  measured  longitudinally  on  duct 
axis  from  initial  station 

airfoil  chord 

mean  inlet  Mach  number 

perpendicular  distance  from  wall  to  point  at  which 
velocity  u  is  measured 

boundary-layer  thickness,  defined  as  perpendicular 
distance  from  wall  to  point  at  which  contribution 
to  integrals  for  6*  and  6  is  negligible 


displacement  thickness 


I  p6 

u  /  u\  \ 

Jo 

vi1  -  >j 

momentum  thickness 


viscosity 

density 

wall  shear-stress  coefficient 

wall  shear-stress  coefficient  obtained  by  equation  of 
Squire  and  Young  (reference  2) 

wall  shear-stress  coefficient  obtained  by  equation  of 
Ludwieg  and  Tillmann  (reference  3) 


turbulent-normal- stress  coefficient 

/  P6  ag 

dissipation  coefficient  II  ^  dy) 

from  experimental  data 
from  correlation  curve 


^  _1_  d(pu' 2) 
)  2q  3x 


h 
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ANALYSIS 

Momentum  Thickness 


In  the  two-dimensional  analysis  of  reference  1,  the  rate  of  growth 
of  the  momentum  thickness  was  calculated  by  the  momentum  equation 


d6  =  JO  _  0  dqA  +  H\ 
dx  2q  q  dx\  2/ 


(1) 


The  Squire  and  Young  (reference  2)  values  of  skin-friction  coefficients 
were  used  with  this  equation  to  calculate  the  rate  of  growth  of  the 
momentum  thickness. 

The  simplest  apparent  procedure  for  extending  equation  (1)  for 
application  to  conical  diffusers  is  to  redevelop  it  in  circular  measure 
by  simple  summation  of  forces  in  the  longitudinal  direction  as  given 
in  reference  h *  for  example.  If  the  assumption  is  made  that  the 
boundary-layer  thickness  is  small  in  relation  to  the  duct  diameter* 
the  following  approximate  equation  results: 


d0  a  22  _  e  da  A  +  H\  £  dD  ,2) 

dx  2q  q  dx  \  2/  D  dx 

This  equation  differs  from  equation  (1)  by  the  presence  of  an  additional 
term  which  provides  for  the  effect  of  the  changing  duct  perimeter.  The 
same  equation  can  be  derived  more  rigorously  by  starting  directly  with 
the  Navier-Stokes  equations. 

In  general*  the  values  of  momentum  thickness  calculated  by 
equation  (1)  for  two-dimensional  flow  and  by  equation  (2)  for  conical- 
diffuser  flow  closely  approximate  the  experimental  values.  Discrep¬ 
ancies*  local  in  character*  between  the  calculated  and  experimental 
gradients  were*  however*  noted.  Figure  1  shows  typical  examples  of 
experimental  and  calculated  momentum-thickness  gradient  and  momentum 
thickness  as  a  function  of  x  as  illustrated  for  the  two-dimensional 
flow  data  from  reference  £.  This  figure  shows  that  over  most  of  the 
range  quite  good  agreement  is  found  between  experimental  and  calculated 
momentum- thickness  gradient  and  momentum  thickness.  Near  the  downstream 
end  of  the  surface*  however*  the  difference  between  experimental  and 
calculated  momentum-thickness  gradient  is  appreciable;  this  difference 
is  responsible  for  the  gradual  divergence  of  the  two  curves  for  momentum 
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thickness.  Several  investigators  (references  6  to  9)  have  observed 
and  commented  upon  discrepancies  of  this  type,  the  occurrence  of  which 
seemed  to  be  associated  with  proximity  to  separation.  Newman  (refer¬ 
ence  10)  suggested  that  this  lack  of  agreement  was  due  to  evaluation  of 
the  momentum  equation  in  terms  of  mean  velocity  and  the  neglect  of  terms 
arising  from  the  presence  of  velocity  fluctuations.  Wallis  (refer¬ 
ence  11)  mentioned  that  the  term  arising  from  the  Reynolds  normal  stress 
is  of  increasing  magnitude  as  separation  is  approached  and  implied  that 
this  term  might  eliminate  these  discrepancies  if  included  in  the  momentum 
equation;  inclusion  of  this  term  gives 


d£  =  IQ  +  -  £  dgA  +  H\ 

dx  2q  2q  q  dx\  2/ 


in  which  the  turbulent-normal-stress  coefficient  T^/2q  is  given  by 
the  equation 


2q  dx 


The  corresponding  momentum  equation  for  conical  diffusers  is 


££  =  12  +  12  _  £a/i  +  H\  _  9 

dx  2q  2q  qdx\  +  2/~D 


An  examination  of  equation  (3)  was  made  by  using  the  two-dimensional 
data  of  reference  5  and  the  results  are  illustrated  in  figure  2.  The 


solid  curve 


T0  TR\ 

—  +  — -  was  obtained  by  solution  of  equation  (3)  in 

2q  2q/ 


which  the  slopes  ^  used  were  obtained  from  a  curve  faired  through 

the  experimentally  measured  momentum  thicknesses.  A  curve  of  the  Squire 
and  Young  skin-friction  coefficient  (T:o/2q)s<Y.  is  also  given  in  this 

figure .  The  decreasing  value  of  the  Squire  and  Young  skin-friction 
coefficient  with  increasing  longitudinal  distance  is  due  to  increasing 
boundary-layer  Reynolds  number.  The  values  of  turbulent-normal-stress 
coefficient  Tjj/2q  evaluated  from  hot-wire  measurements  in  refer¬ 
ence  5  are  included  in  figure  2.  It  is  apparent  that  the  sum  of  the 
turbulent-normal-stress  coefficient  Tp/2q  and  the  skin-friction 
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coefficient  ( tq/2c i)g  y  exceeds  the  value  of  the  sum  obtained  by- 

evaluation  of  equation  (3)  from  the  18-foot  station  to  approximately 
the  22-foot  station. 

Experiments  by  Ludwieg  and  Tillmann  (reference  3)  by  a  hot-plate 
technique  have  indicated  that  the  skin-friction  coefficient  may  diminish 
with  increasing  shape  factor  H,  a  variable  consideration  of  which  was 
beyond  the  scope  of  the  Squire  and  Young  representation  (referenced). 
Values  of  skin-friction  coefficient  ^Tg/2q^  ^  calculated  from  the 

Ludwieg  and  Tillmann  correlation  (reference  3)  are  also  shown  in  figure  2. 
The  sum  of  the  values  of  skin-friction  coefficient  (To/2q)^  ^  and 

the  turbulent-normal- stress  coefficient  Tft/2q  gives  a  good  approxi- 

t0 

mation  to  the  sum  —  +  —  as  evaluated  from  equation  (3)  from  the 

18-foot  to  the  22-foot  station.  The  approximation  ceases  to  be  good, 
however,  beyond  the  22-foot  station,  the  points  departing  greatly  but 
without  apparent  system  from  the  curve.  This  deviation  may  indicate 
significant  incompleteness  of  equation  (3)  in  this  region.  However,  it 
may  be  noted  that  the  ('TrR/^<l)exp  values  also  depart  greatly  and 

apparently  without  system  from  the  trend  that  might  have  been  expected 
on  the  basis  of  preceding  values.  This  departure  suggests  the  alternate 
possibility  that  the  origin  of  the  deviations  may  lie  in  increased 
difficulty  in  determining  the  experimental  value  of  Tj^/2q  at  the 
downstream  end  of  the  surface. 

The  change  in  turbulent-normal-stress  coefficient  with  increasing  x 
is  noted  in  figure  2  to  be  opposite  in  direction  to  that  of  the  skin- 
friction  coefficient  ip  over  this  range,  and  the  sum  of 

(TQ//2q)jj#'p.  and  T^/2q  does  not  differ  greatly  from  values  of 
^To/2q^s  y  for  the  first  several  feet. 

Because  the  pressure-gradient  term  in  equation  (3)  usually  makes  a 
very  substantial  contribution  to  the  rate  of  growth  of  momentum 
thickness,  the  differences  resulting  from  the  choice  of  equation  (l)  in 
which  (To/2q)g  j  is  used  or  equation  (3)  in  which  (T0/ 2q)L,Ta  an<^ 
Tj^2q  are  used  are  not  so  great  as  might  be  surmised  from  examination 
of  figure  2.  A  better  perspective  is  possible  by  reference  to  fig¬ 
ure  1(a),  in  which  the  use  of  equation  (3)  with  (tq /2q)L.T.  anc^  TR/2q 
yields  a  curve  coincident  throughout  most  of  its  length  with  the  experi¬ 
mental  curve  shown  by  the  solid  line  in  comparison  with  the  dashed  line 
which  was  computed  from  equation  (1)  with  (^o/^^S . Y.  •  The  deviation  of 
the  curve  obtained  by  using  equation  (3)  with  (^o/2q)L.T.  anc*  tr/2c1 
is  indicated  at  the  downstream  end  of  the  surface.  Although  the 
numerical  differences  between  the  use  of  (To/2q)s  Y  anc^  ^he  use  ° ^ 
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the  sum  of  (^o/^^L.T.  an<^  tr/^c1  are  indeed  small  in  the  case 

illustrated,  the  concept  of  diminishing  skin-friction  coefficient  with 
increasing  velocity-profile  shape  factor  H  as  indicated  by  (TQ/2q)L  T 

seems  more  logical  than  the  concept  indicated  by  (TQ/2q)g  y  that  the 

skin-friction  coefficient  is  independent  of  the  velocity-profile  shape. 

On  this  basis,  the  assumption  has  been  made  that  the  use  of  equations  (3) 
and  (3)  with  (To/2q)]j>‘j'.  and  TIt/2q  is  more  nearly  correct.  In 
subsequent  reference  to  these  equations  by  number,  the  use  of 
(to/2i)l.T.  an^  tr/2<1  is  intended. 


Boundary-Layer  Shape  Parameter 

The  equation  given  in  reference  1  for  the  rate  of  growth  of  the 
boundary-layer  shape  factor  is 


dH 

dx 


eli.680(H-2.975) 


6  dq  2q 
q  dx  tq 


2.033(H  -  1.286) 


Unlike  the  momentum  equation,  this  expression  is  basically  empiricalj 
moreover,  it  was  developed  for  two-dimensional  application  and  the 
correct  procedure  for  adaptation  to  conical-diffuser  flow  is  obscure. 
Attempts  to  use  equation  (6)  for  the  calculation  of  diffuser  flows  were 
generally  discouraging.  An  example  of  such  a  case  with  two-dimensional 
thickness  parameters  is  shown  in  figure  3.  The  calculated  gradients 
strongly  exceed  the  experimental  gradients  and  the  calculated  values 
of  H  diverge  from  the  experimental  curve  by  a  steadily  increasing 
margin.  This  example  is  typical  of  results  obtained  in  analyses  of  a 
number  of  similar  flows.  In  some  instances,  adjustments  could  be  made 
in  the  numerical  coefficients  of  equation  (6)  which  would  improve  the 
calculation  for  a  given  example  of  diffuser  flow  without  perceptibly 
altering  the  characteristics  of  the  equation  as  applied  to  the  airfoil 
flows  for  which  it  was  originally  developed.  These  cases  were  special, 
however,  and  no  single  set  of  coefficients  was  found  that  offered 
substantial  improvement  in  some  examples  without  detriment  in  others. 

Because  of  the  difficulties  and  ambiguities  in  the  use  of 
equation  (6),  a  comparable  equation  was  sought  through  analysis  of  the 
kinetic -energy  equations  of  the  boundary  layer.  As  a  result  of  this 
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investigation,  the  following  equation,  for  which  the  derivation  for  two- 
dimensional  flow  is  presented  in  appendix  A,  was  obtained: 


dH  =  9  dq  H(3H  -  1)(H  -  1)  T0  ,v  Ts  (3H  -  l)2  m 

dx  ‘  q-dx  2  2q  V  '  2q  2  yu 


In  this  equation  the  term  Ts/2q  is  defined  as 


2q  Jo  2q  dy 


(8) 


Inherent  in  the  derivation  of  equation  (7)  is  the  assumption  that  the 
actual  boundary-layer  velocity  distribution  may  be  approximated  with 
acceptable  accuracy  by. power  profiles.  An  essentially  identical ‘ equation 
in  more  general  form  and  derived  by  a  more  general  process  has  been 
presented  by  Tetervin  and  Lin  (reference  12).  A  comparison  of 
equation  (7)  with  the  comparable  equation  of  reference  12  is  given  in 
appendix  B.  Tetervin  and  Lin  (reference  12)  point  out  that  the  equation 
is  equally  applicable  to  three-  and  two-dimensional  flows. 

The  opportunity  for  determining  whether  equation  (7)  yields  values 
of  dH/dx  in  agreement  with  experiment  when  a  chosen  skin-friction  law 
and  measured  values  of  TS/2q  are  used  is  seriously  restricted  because 
of  the  paucity  of  data  on  the  shear-stress  distribution  in  a  turbulent 
boundary  layer.  In  reference  however,  are  found  hot-wire  data  from  . 
which  the  values  of  TS/2q  may  be  evaluated  for  several  flow  conditions. 
When  the  experimental  values  of  Ts/2q  were  used  together  with 
(To/2q equation  (7)  yielded  shape-factor  gradients  in  good  agree¬ 
ment  with  those  observed  experimentally.  This  agreement  is  illustrated 
in  figure  U,  in  which  the  solid  curve  gives  the  shape-factor  gradient 
derived  from  a  curve  faired  through  the  experimental  data,  and  the  points 
indicated  are  calculated  as  just  described. 

Equation  (7)  as  derived  and  applied  in  the  present  analysis  is  based 
upon  mean  velocities  and  neglects  those  terms  which  would  arise  from 
turbulent  normal  stresses.  The  order  of  agreement  between  -values  from 
equation  (7)  and  from  experiment  has  been  such  as  to  justify  this 
omission  in  calculating  the  growth  of  the  shape  factor  for  all  conditions 
investigated.  Subsequent  to  the  preparation  of  the  current  analysis, 
other  research  workers  at  the  Langley  Laboratory  have  supplied  an 
extended  form  of  the  energy  equation  incorporating  turbulent-normal- 
stress  terms.  Because  both  the  momentum  and  kinetic-energy  equations 
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are  derived  from  the  Prandtl  boundary-layer  equation,  the  inclusion  of 
the  Reynolds  normal-stress  term  in  the  momentum  equation  seems  to 
require  its  use  in  the  kinetic-energy  equation.  It  has  been  found, 
however,  that,  for  the  data  of  reference  5>,  the  net  change  in  calculated 
rate  of  growth  of  shape  factor  arising  from  the  inclusion  of  these 
additional  terms  is  indeed  quite  small  and  in  the  present  investigation 
they  are  omitted.  Appendix  C  gives  the  details  of  this  investigation. 


CORRELATION  OF  STRESS  COEFFICIENTS 


If  equations  (3)  and  (7)  are  to  be  used  for  calculating  boundary- 
layer  development  from  known  values  of  displacement  thickness,  momentum 
thickness,  and  shape  factor  and  from  stream  conditions  outside  the 
boundary  layer,  it  is  necessary  to  express  the  three  stress  coefficients 
to/2c1>  ts/2c1j  and  TR/2q  directly  or  indirectly  in  terms  of  the 

assumed  basic  variables  (nondimensional  pressure  gradient  i  & 

\  q  cbr  . 

boundary-layer  shape  factor  H,  and  boundary-layer  Reynolds  number  Rq] . 
The  fundamental  correctness  of  such  a  representation,  however,  is  ' 
doubtful.  Dryden  (reference  13)  has  pointed  out  that  the  turbulent 
fluctuations  and  turbulent  shear  stress  at  a  point  in  the  boundary  layer 
cannot  be  regarded  as  directly  related  to  the  mean  speed  and  its  deriva¬ 
tives  at  that  point.  The  possibility  exists,  however,  that  such 
representations,  although  fundamentally  inexact,  may  be  adequate  for  the 
purpose  of  calculating  momentum— thickness  and  shape— factor  gradients  from 
equations  (3)  and  (7).  Furthermore, .it  is  hoped  that  a  more  successful 
approximation  for  a  coefficient  might  be  obtained  by  an  integration 
across  the  boundary  layer  than  would  be  possible  by  representing  the 
value  of  the  integrand  at  any  local  point  within  the  boundary  layer.  In 
view  of  this  possibility  of  adequacy  for  the  purposes  intended  and  the 
obvious  utility  of  such  representations  if  found,  an  investigation  to 
find  such  representations  was  undertaken. 

With  regard  to  the  correlations  that  follow  it  is  emphasized  that 
these  correlations  represent  those  which  have  given  the  best  results  in 
the  simplest  form  for  the  data  that  have  been  examined  to  date.  They 
are  strictly  empirical  and  are  regarded  as  interim  representations  only. 
As  more  data  covering  a  broader  range  of  conditions  become  available,  it 
will  be  necessary  to  revise  these  correlations  not  only  for  greater 
numerical  accuracy  but  also  for  improved  correctness  of  representation. 

Skin-friction  coefficient  (TQ/2q)L<T> .-  Ludwieg  and  Tillmann 
(reference  3)  from  their  experiments  with  the  hot-plate  technique 
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reported  their  findings  in  the  form  of  the  following  equation  for  the 
skin-friction  coefficient  in  terms  of  the  Reynolds  number  Rq  and 
shape  factor  H : 


_ 0.123 

\2q/L«T.  ~  Rq 0.268  x  1C)0.678H 


¥ 

f 


This  equation  is  plotted  in  figure  3  for  several  values  of  H.  The 
scope  of  the  experiments  in  reference  3  was  insufficient  to  create  a 
feeling  of  confidence  that  other  variables,  not  considered  in  the 
analysis,  may  not  have  influenced  the  results.  In  view  of  the  good 
results  obtained  in  the  present  investigation  by  using  this  represen¬ 
tation,  however,  and  the  lack  of  definite  evidence  to  the  contrary,  the 
representation  of  skin-friction  coefficient  from  reference  3  has  been 
adopted  for  the  present  analysis. 

Dissipation  coefficient  TS/2q.~  A  representation  of  ms/2q  as  a 

function  of  -  £  H,  and  R0  was  obtained  from  a  detailed  study  of 

data  representing  a  large  variety  of  combinations  of  flow  conditions  and  * 

boundary-layer  parameters.  This  representation  appeared  adequate  for  all 

values  of  0  — .  The  only  directly  measured  values  of  t  /2q  available  9 

dx 

were  those  of  reference  3.  For  other  combinations  of  flow  conditions 
from  references  1,  8,  and  9  and  data  from  conical-diffuser  tests,  xs/2q 
was  obtained  by  solution  of  equation^ ( 7 )  for  experimental  values 

of  _  £  £S  H,  Ro,  and  0  £2.  The  correlation  thus  obtained  is  pre- 
q  dx  °  dx 

sented  in  figure  6. 

Over  the  pressure-gradient  range  up  to  -  ^  ^  =  0.008  the  data  from 

independent  sources  were  sufficient  to  give  reasonable  confidence  in  the 
accuracy  of  the  representation.  The  data  that  were  investigated  at 
higher  pressure  gradients,  although  insufficient  to  establish  the  curves, 
were  not  contradictory  to  the  extrapolation  shown  by  the  dashed  lines  in 
figure  6.  For  study  of  flow  development  at  higher  pressure  gradients, 
these  extrapolations  are  suggested  for  use  until  data  sufficient  to 
establish  their  position  with  greater  accuracy  are  obtained. 

The  precedure  by  which  this  correlation  was  developed  started  with 
extensive  plotting  and  cross-plotting  of  the  coefficient  xs/2q  against 

the  several  basic  parameters  —  —  — ,  H,  and  Ro,  in  this  order,  in 

q  dx 

search  of  trends.  By  assumption  and  trial  of  tentative  algebraic  » 
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representations  of  such  apparent  trends  as  could  be  detected,  a  semblance 
of  order  was  achieved.  A  linear  variation  with  the  nondimensional 
pressure  gradient  was  assumed,  and  an  apparently  systematic  increase 
in  Ts/2q  with  increasing  H  which  appeared  to  be  applicable  over  much 
of  the  pressure-gradient  range  was  found.  This  result  was  cross-checked 
by  giving  an  algebraic  representation  to  the  variation  of  xs/2q  with 
shape  factor  in  order  to  explore  the  variation  with  nondimensional 
pressure  gradient.  With  these  general  relations  blocked  out,  the 
algebraic  representations  were  modified  and  refined  to  incorporate  the 
characteristics  of  low-pressure-gradient  boundary-layer  development  with¬ 
out  substantially  altering  the  relationships  determined  at  larger 
pressure  gradients.  In  these  relations,  the  Reynolds  number  factor  was 
put  in  a  form  parallel  to  that  of  Ludwieg  and  Tillmann  (reference  3) 
for  consistency.  Because  an  orderly  interrelationship  of  the  several 
variables  developed  in  the  process,  it  was  possible  to  effect  appreciable 
improvement  in  the  representation  by  readjustment  of  the  various 
constants.  In  graphing  the  final  equation  the  Reynolds  number  term  and 
zero-pressure-gradient  -  shape-factor  term  were  grouped  with  the  coef¬ 
ficient  Ts/2q  for  clarity  of  representation. 

The  order  of  dispersion  of  the  data  from  the  correlation  finally 
adopted  is  illustrated  in  figure  7.  In  this  figure  the  ratio  of  the 
experimental  value  of  Tg/2q  to  the  correlation-curve  value  is  plotted 
against  the  nondimensional  pressure  gradient  for  five  groups  of  data  dis¬ 
tinguished  by  the  ranges  of  shape-factor  values  into  which  the  data  fall. 

Turbulent-normal-stress  coefficient  Tp/2q.~  Values  of  TR/2q 
calculated  from  hot-wire  data  were  obtained  from  reference  5>  and 
additional  values  for  other  combinations  of  flow  conditions  and  boundary- 
layer  parameters  necessary  for  the  achievement  of  a  correlation  were 
obtained  from  solution  of  equation  (3)  or  (3)  by  use  of  data  from 
references  1,  8,  and  9  and  from  conical-diffuser  tests.  The  correlation 

obtained  is  given  in  figure  8.  The  coordinate  0  22,  although  a  derived 
quantity  and  not  one  of  the  basic  parameters,  is,  however,  a  function  of 
the  basic  variables.  Therefore  ||  is  an  implicit  function  of  the 

basic  variables.  This  parameter  was  adopted  after  a  study  of  plots 
against  a  number  of  variables  and  combinations  thereof.  Although 
suggestions  of  trends  appeared  in  several  of  the  plots,  the  dispersion 
of  the  data,  wide  as  it  is,  was  definitely  less  in  the  plot  finally 
selected.  The  degree  of  certainty  associated  with  the  values  of  the 
various  data  points  varies  considerably  and  depends  upon  the  method  of 
determination  and  the  test  conditions.  Consideration  was  given  to  this 
circumstance  in  constructing  the  curve  given  for  the  coordinates 
selected.  Some  of  the  dispersion  of  the  points  in  figure  8  may  arise 
from  inability  to  fix  the  value  of  TR/2q  to  the  desired  degree  of 
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accuracy  by  the  method  employed.  It  is  equally  probable,  however,  that 
the  dispersion  shown  may  be  indicative  of  the  absence  of  a  significant 
correlating  factor. 


RESULTS  AND  DISCUSSION 


The  practical  application  of  either  equation  (3)  or  (5)  and 
equation  (7)  to  turbulent-boundary-layer  development  must  be  regarded 
as  essentially  empirical.  The  utility  of  the  equations  is  dependent  on 
the  success  achieved  in  attaining  correlations  for  the  factors  ro/2q, 
Ts/2q,  and  TR/2q,  and  the  degree  of  success  is  measured  by  the  accuracy 
obtained  in  applying  the  equations.  In  the  figures  which  follow 
(figs.  9  to  19),  the  accuracy  of  the  analysis,  in  which  the  empirical 
representation  of  the  several  coefficients  is  treated  as  an  integral 
part,  is  demonstrated  by  comparisons  between  experimental  and  calculated 
values  of  momentum  thickness  0  and  shape  factor  H  for  a  represen¬ 
tative  variety  of  examples.  These  demonstrations  are  made  by  two 
procedures.  In  the  first  of  these  procedures,  the  gradients  of  0 
and  H  are  calculated  from  equation  (3)  or  (5)  and  equation  (7)  from 


experimentally  determined  values  of 


0  dq 
q  dx* 


H,  and  R@  at  each  £oint 


of  measurement.  These  gradients  are  plotted  against  distance  along  the 
surface  and  integrated  by  means  of  an  integraph  to  obtain  curves  of  0 
and  H.  These  curves  are  intended  to  show  the  general  nature  of  local 
agreement  between  calculated  and  experimental  values  in  a  manner  as 
free  as  possible  of  errors  arising  from  inaccuracies  in  upstream  regions. 
These  curves  are  identified  by  the  phrase  "local  check."  The  second 
procedure  for  examining  the  accuracy  of  the  method  is  to  calculate  the 
development  of  0  and  H  as  a  function  of  longitudinal  distance  by  a 
step-by-step  integration  in  which  only  the  stream  conditions  at  the 
edge  of  the  boundary  layer  and  the  initial  values  of  0  and  H  are 
used.  Curves  of  0  and  H  as  a  function  of  x  or  X,  obtained  by 
this  process,  are  identified  as  "predicted"  curves  and  are  intended  to 
show  by  comparison  with  the  "local  check"  curves  the  extent  to  which 
cumulative  errors  of  integration  affect  the  result. 


Two-dimensional  flow.-  The  first  of  these  examples  (fig.  9)  is  for 
the  two-dimensional  data  of  reference  £.  The  agreement  obtained  is 
excellent  by  either  method  for  both  0  and  H;  however,  it  should  be 
borne  in  mind  that  these  data,  which  were  much  more  detailed  than  any 
of  the  data  from  other  sources,  were  drawn  upon  heavily  in  formulating 
the  analysis. 


Figure  10  shows  results  as  applied  to  the  data  of  reference  9 
obtained  from  tests  on  a  flat  plate  with  a  boundary  layer  of  the  same 
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general  size  as  that  of  reference  5  but  generated  in  a  channel  differing 
somewhat  in  geometric  detail.  As  the  end  of  the  channel  is  approached, 
the  calculated  values  of  8  deviate  from  the  experimental  values.  Fig¬ 
ure  11  shows  results  as  applied  to  the  very  similar  tests  of  Wieghardt 
(reference  8). 

The  results  applied  to  additional  data  from  reference  9 ,  where  the 
pressure  gradient  on  the  flat  plate  was  adverse  to  only  a  very  slight 
degree,  are  presented  in  figure  12.  Agreement  in  H  is  especially 
good;  the  calculated  values  of  momentum  thickness  0  are  slightly  above 
the  experimental  values  at  the  downstream  end.  In  this  example,  the 
contribution  of  the  pressure  gradient  to  the  rate  of  change  of  shape 
factor  is  very  small,  and  it  is  to  be  expected  that  the  Reynolds 
number  Rq  should  be  the  dominant  influence  in  determining  the  value  of 
shape  factor  H.  That  this  presumption  appears  to  be  valid  is  illus- 
strated  in  figure  13  in  vtfiich  the  experimental  values  of  shape  factor 
from  figure  12  are  plotted  against  the  Reynolds  number  Rq  together  with 
a  curve  obtained  by  a  fairing  of  the  experimental  data  taken  on  a  flat 
plate  at  zero  pressure  gradient,  both  those  of  Ludwieg  and  Tillmann  (ref¬ 
erence  3)  and  also  those  of  Schultz-Grunow  (reference  lit). 

v  As  a  final  case  of  two-dimensional  flow,  figure  Ik  shows  the 

calculation  made  for  the  data  from  the  example  of  boundary-layer  develop¬ 
ment  over  a  two-dimensional  airfoil,  cited  in  reference  1.  Here  again 

*  the  agreement  is  generally  good,  with  the  notable  exception  of  the  last 
point  which  immediately  preceded  separation.  In  this  specific  example, 
the  result  obtained  by  the  present  method  is  unsatisfactory  and  not  so 
good  as  that  shown  for  this  example  in  reference  1. 

Conical-djffuser  flow.-  Four  examples  of  boundary-layer  flow  in  a 
conical  diffuser  are  shown  in  figures  1S>,  16,  18,  and  19.  The  diffuser 
for  these  examples  was  a  12°  cone  with  a  2:1  outlet-to-inlet  area  ratio 
and  a  21-inch-diameter  inlet.  Calculations  were  started  at  the  first 
measuring  station  downstream  of  the  faired  junction  between  the  cylin¬ 
drical  inlet  duct  and  conical  section.  Calculations  upstream  of  this 
point,  although  considered  generally  satisfactory ,  are  not  offered  as 
evidence  of  the  accuracy  of  the  analysis  because  of  the  insufficiency  of 
experimental  data  in  this  region  where  the  pressure  changes  most  rapidly. 
The  first  of  this  group  of  figures  (fig.  l£)  is  for  an  inlet  boundary 
layer  having  a  momentum  thickness  of  0.00217  foot  and  a  shape  factor 
of  1.226  at  a  mean  Mach  number  in  the  approach  duct  of  0.27.  The  second 
(fig.  16)  is  for  the  same  Mach  number  with  an  inlet-boundary-layer 
momentum  thickness  of  0.0112  foot  and  a  shape  factor  of  1.220.  For  both 
of  these  figures  the  agreement  between  theory  and  experiment  is  quite 

*  good  except  at  the  final  points.  The  station  at  which  the  final  points 
were  obtained  actually  was  in  the  upstream  end  of  the  cylindrical 
discharge  duct  rather  than  in  the  diffuser  proper.  Because  of  this* 

*  circumstance,  the  validity  of  using  these  points  to  assess  the  accuracy 
of  the  procedure  is  questionable. 


The  conical-diffuser  flow  illustrated  in  figure  16  is  that  used  in 
figure  3  to  illustrate  the  inadequacy  of  equation  (6).  The  results 
obtained  with  equation  (7)  plotted  to  the  same  coordinates  as  those  of 
figure  3  are  shown  in  figure  17.  The  improvement  in  agreement  as 
obtained  in  equation  (7)  is  evident  from  a  comparison  of  figures  3  and  17 

Although  the  present  analysis  is  for  incompressible  flow,  the 
results  appear  to  be  acceptable  for  moderately  high  speeds,  as  illus¬ 
trated  by  figure  18,  in  which  calculations  are  presented  for  an  initial 
boundary-layer  momentum  thickness  of  0.00ll±8  foot,  a  snape  factor 
of  1.225,  and  a  mean  inlet  Mach  number  of  0.1;3.  The  principal  discrep¬ 
ancy  is  at  the  last  points.  In  figure  19,  calculations  are  shown  for  an 
initial  momentum  thickness  of. 0.0112  foot,  a  sh^pe  factor  of  1.223, 
and  a  mean  inlet  Mach  number  of  0.5U* 


Separation.—  It  is  not  the  purpose  of  this  paper  to  propose  new 

criteria  for  the  calculation  of  the  point  of  separation  of  the  boundary 
layer  nor  to  consider  the  relative  merits  of  such  criteria  as  may  have 
been  proposed  by  others.  It  is  of  interest  to  note,  however,  that  in  the 
process  of  making  the  calculation  of  the  boundary-layer  development  for 
the  data  of  reference  5  the  trend  of  the  calculated  skin-friction- 
coefficient  curve  was  observed  t a  be  toward  zero  at  a  point  at  which 
separation  had  been  established  by  direct  measurements  in  the  original 
experiment.  This  trend  is  illustrated  in  figure  20.  Although  the 
Ludwieg  and  Tillmann  correlation  (reference  3)  makes  the  skin-friction 
coefficient  (^o/2q)L.T.  asymptotic  to  zero  as  the  shape  factor  H 

becomes  infinite,  at  the  last  point  at  which  measurements  were  made 
prior  to  the  occurrence  of  separated  flow,  it  appears  that  the  rate  of 
increase  of  the  shape  factor  H  with  distance  is  itself  increasing  at 
a  rate  sufficient  to  give  the  curve  of  skin-friction  coefficient 
(To/2q)k  ip  a  trend  suggestive  of  a  finite  intercept  on  the  X-axis. 

APPLICATION 


The  preceding  discussion  has  been  concerned  with  the  accuracy  of  the 
analysis  which  was  examined  by  application  to  a  representative  group  of 
examples  by  two  methods.  Inasmuch  as  the  principal  anticipated  utility 
of  the  current  analysis  would  be  the  prediction  of  the  boundary-layer 
development,  the  accuracy  of  the  analysis.,  as  demonstrated  by  the  second 
procedure  (performing  a  stepwise  integration  using  only  the  initial 
values  of  0  and  H  and  stream  conditions  outside  the  boundary  layer), 
is  of  particular  interest.  The  necessity  of  repeated  solution  of 
either  equation  (3)  or  (5)  and  equation  (7),  however,  renders  this  step¬ 
wise  calculation  especially  inconvenient.  In  order  to  facilitate 


NACA  TN  2U78 


* 


w 


% 


1$ 


execution  of  this  procedure,  the  terms  of  either  equation  (3)  or  (5>) 
and  equation  (7)  have  been  grouped  and  graphed  in  a  convenient  manner. 

The  analysis  as  so  organized  is  given  by  the  following  equations  together 
with  associated  figures: 

e  =  1  +  11  (10) 

dx 

For  two-dimensional  flow: 


=  III  +  IV  +  V  (11) 

For  conical  diffuser  flow: 


~=III  +  IV+V-|^  (12) 

dx  D  dx  v  ' 

Algebraic  expressions  corresponding  to  the  graphs  of  terms  I,  II,  III, 
IV,  and  V  are  given  in  appendix  D. 


The  terms  I  and  II  in  equation  (10)  are  obtained  from  figures  21 
and  22,  respectively,  as  functions  of  the  nondimensional  pressure 

gradient  -  “  ^  and  the  shape  factor  H.  Term  I  is  read  directly  from 

figure  21.  Term  II  is  obtained  by  dividing  the  ordinate  of  figure  22 
by  Rg 0*268 .  Although  the  curves  of  figure  22(a)  appear  to  have 

9  dq 

a  common  intercept  at  —  ^  =  0,  some  displacement  of  the  intercept 

actually  occurs  for  different  values  of  H,  and  this  displacement  is  of 
importance  when  developments  at  very  small  pressure  gradients  are  to  be 
calculated.  An  enlargement  of  the  graph  for  small  pressure  gradients 
is  given  in  figure  22(b),  and  the  zero-pressure-gradient  intercepts  are 
plotted  as  a  function  of  shape  parameter  in  figure  22(c).  Term  III  for 
equations  (11)  and  (12)  is  given  in  figure  23  as  a  function  of  the 

nondimensional  pressure  gradient  -  ^  ^  and  shape  factor  H.  Term  IV 

is  (^o/^L.T..  an^  ■'■s  represented  in  figure  5  as  a  function  of  the 

Reynolds  number  Re  and  shape  factor  H.-  Term  V  is  xft/2q  and  is 
represented  in  figure  8  as  a  function  of  the  shape-factor-growth  coef- 
dH 

ficient  9  ^  which  is  made  available  by  the  solution  of  equation  (10). 

For  calculating  the  boundary-layer  growth  in  a  conical  diffuser  the  duct 
dimensions  must  be  known  and  equation  (12)  should  be  used. 
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CONCLUSIONS 


A  procedure  based  on  the  kinetic-energy  equation  and  an  extended ‘ 
form  of  the  momentum  equation  has  been  devised  for  calculating  the 
development  of  turbulent  boundary  layers,  in  adverse  pressure  gradients, 
for  that  class  of  flows  for  which  the  fluid  density  at  all  points  and 
total  pressure  outside  the  boundary  layer  are  invariant.  In  the  develop¬ 
ment  of  this  procedure  an  effort  was  made  not  only  to  arrive  at  an  ana¬ 
lytical  form  that  would  allow  examination  of  the  significance  of  the 
physical  quantities  involved  but  also  to  achieve  a  high  degree  of  con¬ 
sistency  with  the  more  recent  results  of  .turbulent-boundary-layer 
research. 

Correlations,  which  are  essential  to  the  execution  of  the  method, 
of  several  elements  of  the  basic  equations  in  terms  of  quantities 
regarded  as  controlling  parameters  are  presented.  These  correlations 
are  strictly  interim  empirical  relationships  both  as  to  numerical  values 
and  as  to  the  variables  involved.  It  is  to  be  expected  that,  with  the 
acquisition  of  more  data  covering  a  broader  range  of  conditions,  improve¬ 
ments  will  be  made  in  both  the  functional  nature  and  the  accuracy  of  the 
correlations.  A  comparable  and  associated  refinement  of  the  basic 
equations  is  likewise  anticipated  as  more  data  are  obtained. 

Predictions,  by  this  method,  of  turbulent-boundary-layer  development 
are  compared  with  experimental  results  from  several  sources  ,f or  a  number 
of  cases  of  flow  over  flat  plates  and  airfoils  and  in  conical  diffusers. 
In  the  range  of  boundary-layer  flow  short  of  separation,  the  agreement 
of  calculated  values  with  experimental  values  was,  in  most  cases,  quite 
satisfactory;  in  some  instances,  however,  definite  disagreement  between 
the  calculated  and  experimental  results  was  noted.  It  is  believed, 
however,  that  good  agreement  has  been  obtained  in  enough  instances  to 
justify  continuation  of  effort  along  the  present  lines,  particularly 
with  respect  to  improvement  of  the  correlations  and  refinement  of  the 
equations. 


Langley  Aeronautical  Laboratory 
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APPENDIX  A 

DERIVATION  OF  EQUATION  FOR  RATE  OF  GROWTH  OF  BOUNDARY-LAYER 
SHAPE  FACTOR  (EQUATION  (7)) 


The  two-dimensional  equation  for  the  loss  of  kinetic  energy  per 
unit  length,  normal  to  the  direction  of  flow,  is  given  in  reference  7  as 


The  kinetic-energy  thickness  of  the  boundary  layer  is  defined  as 


If  equation  (Al)  is  divided  by  U  and  the  definition  for  the 
kinetic-energy  thickness  63  is  substituted  in  the  right-hand  side  of 
the  equation,  the  following  equation  results: 


-nU 

^  .  1 
t  —  dy  =  — 

ay  u3 


du3§3 

dx 


q 

By  substitution  of  the  definition  for  Ts/2q 


,(A2) 


•  is.  r6.^, 

2q  Jo  2q  dy  7 

equation  (A2)  becomes 

ls=li  dU^63 

2q  2  u3  dx 


(A3) 
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By  transposing  terms  and  substituting  the  representation  given 
for  63,  the  following  equation  is  obtained  from  equation  (Alt): 


d63  d  /  UH  0  ts  3  S3  dq 
"33T  "  dx  ^  3H  -  lj  2q  ~  2  q  dx 


Differentiating  the  left-hand  side  of  this  equation  and  transposing 
terms  yields  an  equation  of  the  following  form: 


6 


dH 

dx 


=  H(3H  -  1) 


d0 

dx 


3H  -  l(  _  3  £3  dq\ 
liH  y  2q  2  q  dxy 


(A6) 


Substituting  the  momentum  equation  (equation  (l))  and  the  definition  for 
63  (equation  (A5))  into  this  relation  yields  the  kinetic-energy  equation 
for  the  rate  of  growth  of  boundary-layer  shape  factor  H  in  its  final 
form  for  the  present  analysis 


dH  _  e  dq  H(H  -  1)(3H  -  1)  T0  ,  ,  Ts  (3H  -  l)2  .  < 

_____ - - -  +  _  H(3H  -  1)  -  ^  ^ -  (A7) 
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APPENDIX  B 


COMPARISON  OF  EQUATION  (7)  OF  THE  PRESENT  ANALYSIS 
WITH  EQUATION  (13)  OF  REFERENCE  12 


Equation  (13)  in  reference  12  is  given  as 


dH  £  dU 
H  dx  '  "  "  U  dx 


H(H  -  1)(3H  -  1)  +  (3H  -  l)|fi  +  —  2"~  Jq  *P  ||  qj  — 


VQ  (H  ±  1)(3H  -  1 

u  IT 


with  the  following  accompanying  definitions: 


t  _  y  _  p 
1  "  F  -  vTJ 


g  =  7- 

To 


Vq  is  the  velocity  perpendicular  to  the  Wall  at  the  wall  (y  ■  0), 
and  p  is  an  exponent  related  to  the  profile  shape.  Substitution 
of  these  definitions  in  equation  (Bl)  and  use  of  the  relation 


0  dU  =  _9_  da 
U  dx  2q  dx 
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gives  the  following  equation: 


0  dH  ,  .  e  dg  H(3H-1)(K-1)  *  22  „(3H  -  1) 
dx  q  dx  2  2q 


To  (3H  -  i)2  P1  u  _0  7  +  12  (H ■+  1)(3H  - 


6  U 


The  relation 


u  I2-  dy  = 


T^dy 

ay 


may  be  nondimensionalized  in  the  following  manner: 


f1  u  3,  .  y  _ 

o  DsS  6~ 


to  y/6  6 


Substitution  of  this  relation  in  equation  (B2)  yields  an  equation  of  the 
following  form: 


dH  _  9  dq  H(3H  ~  1)(H  -  1)  ,  ^0  ,_H  l} 
dx  q  dx  2  2q 

T0  (3H  -  l)2  f1  t  ^  y  v0  (H  +  1)(3H 

2i  2  Jo  x°d  5  ■  »  r" 


Substituting  the  definition  for  Ts/2q 


o  2M  6 

o 


■a. 
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f 

into  equation  (B3)  and  noting  that  the  term  Vq/U  becomes  zero  for  the  , 

solid-wall  case  being  considered  herein  yields  I 


0®=J  dq.H(3H  -  1)(H  -  1)  +  10  H(3H  _  1} 
dx  q  dx  2  2q 


Is  (3H  -  l)2 
2q  2 


The  result  obtained  is  the  same  as  equation  (7)  of  the  present  analysis. 
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APPENDIX  C 

EXTENDED  ENERGY  EQUATION 


The  energy  equation  as  extended  to  include  turbulent-nonnal-stress- 
coefficient  terms  is: 


dH  =  _  6  dq  H(H  -  1)(3H  -  1)  +  ^0  R(  _  ^  (3H  -  l)2 

dx  q  dx  2  2q  k  '  2q  2 


IS 

2q 


H(3H  -  1) 


in  which,  as  before. 


TR  i  p6  a(Pu'2)  , 

“Es r * 

and  the  new  term  Tjj»/2q  is  defined  as 


It  may  be  noted  that  equation  (Cl)  differs  from  equation  (7)  by  the 
addition  of  the  terms  in  Td/ 2q  and  Tj^*/2q.  The  sum  of  these  two 
terms,  divided  by  6,  therefore,  constitutes  the  correction  to  the  calcu¬ 
lated  dhape— factor  gradient  which  is  to  be  made  to  provide  for  turbulent 
normal  stresses.  Figure  2k  shows  the  experimental  data  of  reference  5 
and  the  shape-factor  development  calculated  from  equation  (7),  trans¬ 
ferred  from  figure  9.  (The  curve  of  dH/dx,  which  was  integrated  to 
obtain  the  curve  of  H  in  figure  9>  has  been  added  for  more  complete 
comparison.)  The  dashed  curve  of  figure  2lj(a)  has  been  obtained  by 
adding  the  turbulent-normal-stress  corrections  to  the  solid  curve,  and 
the  integration  of  this  dashed  curve  constitutes  the  dashed  curve  of  H 
against  x  in  figure  2U(b).  In  general,  the  departures,  either  in  rate 
of  change  or  value  of  the  shape  factor,  are  not  great.  In  interpreting 


2h 
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figure  21*,  the  significant  item  is  the  deviation  of  the  dashed  curves 
from  the  solid  curves,  not  the  general  position  of  either  relative  to 
the  data,  because  the  empirical  representation  of  TS/2q  was  established 
without  regard  to  the  turbulent— normal— stress  terms. 
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matical  fitting  of  the  correlations  obtained  and  should  not  under  any  circumstances  be  regarded 
as  having  any  other  significance  or  applicability  outside  the  range  of  the  curves  shown. 
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(h)  Momentum  thickness. 

Figure  1.-  Comparison  of  calculated  and  experimental  momentum- thickness 

gradient  and  momentum  thickness. 
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Figure  2.-  Variation  with  longitudinal  distance  of  skin- friction 
coefficient  and  turbulent- normal-stress  coefficient. 


Experimental 


Calculated  (reference  1) 
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Figure  4.-  Comparison  of  calculated  and  experimental  shape-factor 
gradients.  Calculated  values  obtained  from  equation  (7)  with 

(4.i. 


003  , 


Reynolds  number  Rg  and  shape  factor 
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A  Reference  1 
O Reference  5 
^Reference  8 

□  Reference  9 
A  Diffusers  M  =  0.27 
(A  Diffusers  M  =  0.43  and  M  ■  0.54 


gradient 


ob  ! 


NACA  TN  21*78 


(a)  Local  check. 

Figure  9.-  Variation  with  longitudinal  distance  of  experimental  and 
calculated  values  of  momentum  thickness  and  shape  factor. 
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(a)  Local  check. 

Figure  10.-  Variation  with  longitudinal  distance  of  experimental  and 
calculated  values  of  momentum  thickness  and  shape  factor. 
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X,  ft 

(b)  Predicted. 


Figure  10.-  Concluded. 


ft 


(a)  Local  check. 

Figure  11.-  Variation  with  longitudinal  distance  of  experimental  and 
calculated  values  of  momentum  thickness  and  shape  factor. 
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Figure  12.-  Concluded. 


(b)  Predicted. 
Figure  14.-  Concluded 


□  Experimental 
—  Calculated 
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(t)  Predicted. 
Figure  l6.-  Concluded 


Figure  17* -  Comparison  of  calculated  (present  method^  equation  (7))  and 
experimental  shape-factor  gradient  and  shape  factor  for  conical- 
diffuser  flov  for  Mach  number  0.27. 
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(t>)  Predicted. 
Figure  18.-  Concluded 
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Figure  19. -  Concluded. 
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Figure  20.-  Variation  of  skin- friction  coefficient  (— \  with  longi 

V  */  L.  T  • 

tudinal  distance  for  two-dimensional  data,  of  reference  5- 
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Figure  22.-  Term  II  x  RqQ-268  of  equation  (10)  expressed  as  a  function 

of  the  nondimensional  pressure  gradient  -  -  ^  for  constant  values 

q  clx 

of  shape  factor  H. 
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